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ABSTRACT 

We explore the infrared M band (4.7 pm) spectrum of the class I protostar L1489 IRS in the Taurus 
Molecular Cloud. This is the highest resolution wide coverage spectrum at this wavelength of a low mass 
protostar observed to date (R =25,000; Av =12 km s" 1 ). A large number of narrow absorption lmes 
of gas phase 12 CO, 13 CO, and C 18 0 are detected, as well as a prominent band of solid CO. The gas 
phase 12 CO lines have red shifted absorption wings (up to 100 km s *), which likely originate from warm 
disk material falling toward the central object. Both the isotopes and the extent of the CO line wings 
are successfully fitted with a contracting disk model of this evolutionary transitional object (Hogerheijde 
2001). This shows that the inward motions seen in millimeter wave emission lines continue to within 
~0.1 AU from the star. The amount of high velocity infalling gas is however overestimated by this 
model, suggesting that only part of the disk is infalling, e.g. a hot surface layer or hot gas in magnetic 
field tubes. The colder parts of the disk are traced by the prominent CO ice band. The band profile 
results from CO in ’polar’ ices (CO mixed with H 2 0), and CO in ’apolar’ ices. At the high spectral 
resolution, the ’apolar’ component is, for the first time, resolved into two distinct components, likely 
due to pure CO and CO mixed with C0 2 , 0 2 and/or N 2 . The ices have probably experienced thermal 
processing in the upper disk layer traced by our pencil absorption beam: much of the volatile apolar 
ices has evaporated, the depletion factor of CO onto grains is remarkably low (~7%), and the C0 2 
traced in the CO band profile was possibly formed energetically. This study shows that high spectral 
resolution 4.7 fi m observations provide important and unique information on the dynamics and structure 
of protostellar disks and the origin and evolution of ices in these disks. 

Subject headings: dust, extinction — Infrared: ISM — ISM: molecules stars: formation stars: 
individual (L1489 IRS)— planetary systems: protoplanetary disks 


1. INTRODUCTION 

In the process of low mass star formation, a mixture of 
gas, dust, and ices accumulates in protostellar envelopes 
and disks. The fate of this molecular material is diverse. 
Most of it will fall toward the protostar and dissociates in 
the inner disk region or stellar photosphere. Some ma- 
terial will be blown away and destroyed by the stellar 
wind. Some may also survive and be the building mate- 
rial for comets and planets. Major aspects of this compli- 
cated process are not well understood, and poorly obser- 
vationally constrained. For example, do the ices that form 
comets still resemble ices of the original pristine molecu- 
lar clouds or are new ices of different composition being 
formed in the envelope or disk? The type of ices being 
formed depends on the composition of the gas that ac- 
cretes onto grains. Reducing environments produce H 2 0- 
rich (‘polar’) ices, while in cold inert environments ‘apolar’ 
ices rich in CO, N 2 and 0 2 can be formed (Tielens & Hagen 
1982). Depending on the composition, ices evaporate be- 
tween temperatures of 18 and 90 K. Also, heat can change 
the solid state structure of ices by for example crystalliza- 
tion. Energetic particles (e.g. cosmic rays) and ultraviolet 


(UV) radiation are able to initiate reactions in ices and 
form new species. Dynamics and shocks within disks may 
be able to destroy ices as well. 

Clearly, to determine the relative importance of these 
ice formation and destruction processes, knowledge of the 
physical conditions and structure of envelopes and disks is 
crucial. Much theoretical and observational work on this 
topic has been done over the last ^10 years. Molecular gas 
was detected in a suite of protostellar disks by millimeter 
wave observations sensitive to emission over radii of sev- 
eral hundred AU (Dutrey, Guilloteau, & Guelin 1997; Thi 
et al. 2001). Gas phase abundances were found to be 
reduced by factors of 5 to several 100, depending on the 
source and the sublimation temperature of the molecules. 
Models of disk mid-planes indeed show high depletions be- 
cause of the formation of icy mantles on grains (Aikawa 
et al. 1997; Willacy et al. 1998). The predicted deple- 
tions were in fact higher than observed and thus desorption 
mechanisms are needed to explain the millimeter wave ob- 
servations (Goldsmith, Langer, & Velusamy 1999). It was 
realized that the outer parts of disks are heated more effi- 
ciently when they are flared (Kenyon & Hartmann 1987). 
Thus, by the influence of the stellar radiation a layer with 
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‘super-heated’ dust is formed in which molecules have been 
dissociated (Chiang & Goldreich 1997). The layer below 
that is warm enough to evaporate the ices, but not dis- 
sociate the released molecules. The importance of this 
warm layer, and relative gas phase molecular abundances, 
depends strongly on how effective ice desorption mecha- 
nisms are (Willacy et al. 1998). Recent studies indicate 
that desorption by UV and X-ray photons may be strong 
enough to explain observations of molecular gas by mil- 
limeter wave telescopes (Willacy & Langer 2000; Najita, 
Bergin, & Ullom 2001). This idea is confirmed by multi 
transition molecular line observations which indicate tem- 
peratures (> 20 - 40 K) and densities that typically occur 
in this warm layer (van Zadelhoff et al. 2001). In a 
third layer, the disk mid-plane, cold, dense conditions pre- 
vail, resulting in extreme depletions of gas phase molecules 
many orders of magnitude larger than in quiescent dense 
molecular clouds. Indeed, recent absorption line obser- 
vations failed to detect gas phase CO in the edge-on disk 
around the protostar Elias 18, thus indicating an enormous 
depletion in the mid-plane (Shuping et al. 2001). 

In this Paper we report high spectral resolution 
(R =25, (WO) 4.7 fim M band observations of the obscured 
protostar L1489 IRS (IRAS 04016+2610) in the Taurus 
Molecular Cloud. L1489 IRS is a low luminosity object 
(3.7 L 0 ), with a spectral energy distribution resembling 
that of an embedded class I protostar (Kenyon, Calvet, 
& Hartmann 1987). Detailed millimeter wave line and 
continuum studies show that L1489 IRS is surrounded by 
a large, 2000 AU radius rotating thick disk-like struc- 
ture (Hogerheijde et al. 1998; Saito et al. 2001), 
rather than an inside-out collapsing envelope (Hogerhei- 
jde & Sandell 20W). The rotation is sub-Keplerian, and 
the disk as a whole is contracting. Thus it was suggested 
that L1489 IRS represents a short-lasting (2xl0 4 yr) tran- 
sitional phase between embedded YSOs that have large 
envelopes and small (few hundred AUJ rotationaUy sup- 
ported disks, and T Tauri stars which have no envelopes 
and fully rotationaUy supported 500-800 AU size disks 
(Hogerheijde 2001). This circumstellar (or circumbinary: 
Wood et al. 2W1) disk is seen close to edge-on (60 to 
<90°) in scattered light images (Whitney et al. 1997; 
Padgett et al. 1999). A CO outflow emanates from the 
object (Myers et al. 1988) with Herbig Haro objects lying 
along it (G6mez, Whitney, & Kenyon 1997). 

Low spectral resolution infrared observations of L1489 
IRS show that deep H 2 O (Sato et al. 1990) and CO 
(Chiar et al. 1998; Teixeira, Emerson, &; Palumbo 1998) 
ice bands are present along the line of sight. In this Paper 
we will use the newly available spectrometer NIRSPEC at 
the Keck II telescope to obtain high resolution M band 
spectra (Av =12 km s _1 at 4.7 /zm) of this source. The 
large array of NIRSPEC allows both the vibrational band 
of solid CO and the surrounding ro- vibrational transitions 
of gas phase CO to be observed in the same high reso- 
lution spectrum. This offers a new view on this system, 
both on the origin and evolution of ices and the interrela- 
tionship of gas and ices as well as on the kinematics and 
structure of the young, contracting, close to edge-on disk. 
It is a unique view, because infrared absorption line stud- 
ies trace all gas and solid state material at all radii from 


the star, while present day millimeter wave observations 
are limited by their relatively low spatial resolution (>100 
AU). Thus one of the questions that will be answered in 
this Paper is whether the large scale inward motions seen 
in millimeter wave emission lines continue to smaller radii 
(1 AU or less) from the star. 

Previous studies have already shown that rich astro- 
physical information can be obtained from high spec- 
tral resolution observations in the atmospheric M band. 
Mainly massive, luminous protostars were observed, how- 
ever (Mitchell et al. 1990) , and the few observations of low 
mass protostars cover a small wavelength range contain- 
ing only a few gas phase lines and not the solid CO band 
(Shuping et al. 2001; Carr, Mathieu, & Najita 2001). 
This is Paper I in a series on high resolution M band spec- 
troscopy of protostars, initiated by the availability of the 
NIRSPEC spectrometer at Keck II with which weak, low 
mass protostars can be routinely observed at high spectral 
resolution over a large wavelength range covering both the 
solid and gas phase CO features. 

The reduction of the long slit spectra is discussed in §2. 
In §3.1, we analyze the observed gas lines, which at this 
resolution even give dynamical information. We use stan- 
dard curve of growth and rotation diagram techniques to 
get a first idea of gas column and temperatures. In order 
to analyze the solid CO band profile in this line of sight, 
a detailed discussion of available laboratory experiments 
of solid CO is given in §3.2. In §4.1 we apply the infalling 
disk model of Hogerheijde (2001) to explain the observed 
gas phase 12 CO and 13 CO line profiles, and constrain the 
physical conditions and structure of the disk. The possi- 
bility of binarity is briefly discussed in §4.2. The gas phase 
analysis is linked to the the solid CO results to determine 
the origin and thermal history of solid CO in §4.3. We 
conclude with suggestions for future work in §5. 

2. OBSERVATIONS 

The infrared source L1489 IRS was observed with the 
NIRSPEC spectrometer (McLean et al. 1998) at the Keck 
II telescope atop Mauna Kea on the nights of 28/29 and 
29/30 January 2001. The sky was constantly clear and dry, 
and the seeing was reasonable (~ 0.5 — 0.7" at 2.2 fim). 
NIRSPEC was used in the echelle mode with the 0.43 x 24" 
slit, providing a resolving power of R = A / AX = 25, 000 
(~12 km s“ 2 ) 6 with three Nyquist sampled settings cov- 
ering the wavelength range 4.615-4.819 fim in the atmo- 
spheric M transmission band. 

The data were reduced in a standard way, using IDL 
routines. The thermal background emission was removed 
by differencing the nodding pair. Each nodding position 
was integrated on for 1 minute, before pointing the tele- 
scope to the other nodding position. The positioning had 
to be done by hand because condensations on the dewar 
window prevented us from using the image rotator, and 
the instrument had to be used in the non-standard ‘sta- 
tionary guiding mode’. In this mode, the condensations 
were well removed by subtracting nodding pairs, although 
the sensitivity of Keck/NIRSPEC was reduced by several 
magnitudes. 

The most critical step in the reduction of this data is 
correction for atmospheric absorption features, and sepa- 


This nominal spectral resolution was verified by measuring the Gaussian width of absorption lines in a variety of astrophysical sources. 
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Fig 1.— Observed, unsmoothed, R=25,000 M band spectrum of L1489 IRS corrected for object and earth velocity (43 km s _1 ) on optical 
depth scale, with gas phase line identifications of 12 CO (deepest lines) and 13 CO ('+’ ^mbols above sp^trum) determined from the HI- 
TRAN database (Rothman et al. 1992). The diamonds indicate the expected position of C ls O absorption lines, with the R(4) and R(5) line 
detections labeled. The prominent broad absorption between 2123-2149 cm" 1 can be fully ascribed to solid CO. Wavelength regions with 
poor atmospheric transmission are not plotted. 


rate telluric and interstellar CO absorption features. For 
our observations of LI 489 IRS this proved to be relatively 
easy, because at the time of our observations the deep in- 
terstellar CO lines are shifted by as much as 43 km s~* 
with respect to telluric CO lines. The low airmass of the 
source, 1.00-1.05, also made this step easier. The standard 
stars HR 1380 (A7V) and HR 1497 (B3V) were observed 
at similar airmass, and their spectral shape and hydrogen 
absorption features were divided out with Kurucz model 
atmospheres. A good telluric correction was achieved, al- 
though residuals are seen near deep atmospheric lines. To 
be safe, we therefore removed the parts of the spectrum 
which have less than 50% of the maximum transmission in 
each setting. This does not affect most of the features in 
the L1489 IRS spectrum, because the velocity shift of 43 
km s” 1 fully separates them from telluric features. The 
final signal-to-noise value on the unsmoothed spectrum 
is ~ 45, after integration times of 11, 14, and 20 min- 
utes on each of the three settings of this M = 4.8 mag- 
nitude (Kenyon et al. 1993) source. The different set- 
tings were wavelength calibrated on the atmospheric CO 
emission lines, and subsequently the three settings were 
combined by applying relative multiplication factors. We 
have not attempted to flux calibrate the spectrum, since 
we are interested in absorption features only. 


3. RESULTS 

The fully reduced echelle spectrum of L1489 IRS shows, 
in great detail, many deep narrow absorption lines of gas 
phase 12 CO and 13 CO, and a few weak lines of C 18 0 
(Fig. 1). These lines were identified, using the line frequen- 
cies in the HITRAN catalogue (Rothman et al. 1992). 
The broad absorption feature between 2122-2149 cm -1 can 
be attributed to the stretching vibration mode of 12 CO in 
ices along this line of sight. 

In order to analyze these gas and solid state absorption 
features, a shallow, second order polynomial continuum 
was applied to derive the optical depth spectrum. The 
solid CO band was then analyzed using available labora- 
tory experiments (§3.2). The derivation of physical param- 
eters from the gas phase lines is highly model dependent. 
First, we will derive temperatures and column densities 
using the standard curve of growth and rotation diagram 
techniques (§3.1). Then, we will independently test an 
astrophysically relevant power law model in §4.1. This in- 
formation is combined in §4.3 to discuss the origin and 
thermal history of the solid CO seen in this line of sight. 

3.1. Gas Phase CO 

The 12 CO lines have a complicated profile. Deep lines 
are present at a velocity of -1-43 km s” 1 with respect to 
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earth, which is the systemic velocity of L1489 IRS at the 
date these observations were done (taking i^ 8r = 5 km s~ l 
from millimeter emission lines; e.g. Hogerheijde & Sandell 
2000). Each of the 12 CO lines is accompanied by a fairly 
prominent wing on the red shifted side (Fig. 2). As an 
exploring step in the analysis, we decomposed the main 

12 CO component and its wing by fitting two Gaussians. 
They are separated by on average 23±6 km s” 1 , and the 
main feature and its wing have widths of FWHM=20±3 
km s" 1 and 53±17 km s -1 respectively. The main 12 CO 
feature at the systemic velocity is resolved. Excess absorp- 
tion is visible on the blue and red sides with respect to 

13 CO (Fig. 2). The 13 CO lines have Gaussian shapes with 
a width equal to the instrumental resolution (FWHM==12 
km s~ ! ). The main absorption feature of the 12 CO lines 
is therefore only in part responsible for the same gas seen 
in 13 CO. 



Fig. 2. — Comparison of the 12 CO and 13 CO (thick gray line) 
spectra of L1489 IRS. The spectra shown are averaged over the ob- 
served 12 CO P(6)-P(15) lines and all 13 CO lines in order to increase 
the signal-to-noise. They are corrected for the source and earth ve- 
locity (43 km s -1 ) and the 13 CO spectrum was multiplied by 3,64 
to facilitate comparison. 

We derived equivalent widths for the isotopes and 12 CO 
components (Table 1), and applied a standard curve of 
growth technique to calculate column densities for each J 
level. The comparison of equal J levels of C 18 0, 13 CO, 
and main 12 CO component provides a handle on the in- 
trinsic line width bo (=FWHM/2vln2), in the assumption 
that all material absorbs at the same velocity (however, 
see §4.1). One also has to assume that the isotope ratios 
are constant along the line of sight ( 12 CO/ 13 CO=80 and 
12 CO/C 18 O=560; Wilson & Rood 1994). The equivalent 
widths of C 18 0 and 13 CO are then simultaneously fit at 
&d> 0.8 km s” 1 . Lower bo significantly (> 3 a) underesti- 
mates the 13 CO lines, with respect to C 18 0. As mentioned 
above, the main 12 CO component is clearly contaminated 
by gas not seen in 13 CO, both on the blue and red shifted 
sides. As a first order correction we lowered the 12 CO 
equivalent widths in Table 1 with 40%, which corresponds 
to unresolved lines at the observed peak optical depth. 
Then, we find that all isotopes are best fit simultaneously 
at &d=1*4±0.1 km s” x . This would be an upper limit if 


the correction for contamination of the 12 CO lines were 
too small. Hence, this curve of growth analysis of the CO 
isotopes shows that, in the assumption that all gas absorbs 
at the same velocity, bo is limited to 0.8< &d< 1 *5 km s” 1 . 

With the column densities per J level at hand, a rota- 
tion diagram was constructed for the 13 CO lines (Fig. 3) to 
derive the total column density and temperature at a num- 
ber of allowed bo values (Table 2). Clearly, the rotation 
diagram shows a double temperature structure, much re- 
sembling that of high mass objects (Mitchell et al. 1990): 
cold (T ~ 15 K), and warm gas ( T ~ 250 K) are present 
along the same line of sight. The column density of the 
cold component toward L1489 IRS is a particularly strong 
function of bo , increasing by an order of magnitude from 
0.7 to 1.3 km s" 1 . An independent CO column of 1.4xl0 19 
cm -2 can be estimated from Ay =29 (Myers et al. 1987), 
assuming the dust and gas are co-spatial . This would sug- 
gest bo is in the 1.0-1. 3 km s~ x range. 



Fig. 3. — Rotation diagram of I3 CO, with the column densities on 
the vertical axis calculated from the curve of growth with &d= 1*0 
km s” 1 (dots: P branch; closed triangles: R branch; open triangles: 
R branch upper limits). The solid line corresponds to tempera- 
tures Trot— 16 K, and 250 K (Table 2) in the Boltzmann equation. 
Open diamonds connected with a dotted line represent the rotation 
diagram for our disk model, which is different from the curve of 
growth analysis because of a lower assumed &d (0.1 km s _1 ) and 
the presence of a velocity gradient (§4.1), although both models fit 
the observed 13 CO lines well. 

For the highly red shifted 12 CO wing a rotation dia- 
gram is constructed as well. Here, we assume that the 
absorption is optically thin, because the wings are absent 
in X3 CO, as evidenced by the high signal-to-noise average 
line profile (Fig. 2). The rotational temperature of this 
highly red shifted gas is similar to that of the warm 13 CO 
component (250^^ K), however the column is an order 
of magnitude lower (Table 2). Finally, assuming that the 
absorption in the wings of the resolved 12 CO main com- 
ponent is optically thin, the column of this blue and red 
shifted low velocity gas (within 10 km s _1 of the systemic 
velocity) is about 25% of the column of the high velocity 
red shifted gas at a similar, although poorly determined, 
temperature. 

Summarizing, this basic analysis indicates that, from 
13 CO, both a large column of cold 15 K gas and a sig- 
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Table 1 

Equivalent Widths 


transition 

12 co 

10 -3 cm 1 
main wing 

transition 

13 co 

W v 

10 -3 cm -1 

transition 

c l8 o 

W u 

10“ 3 cm" 1 

R(5) 

88 (11) 

70 (22) 

R(17) 

< 3 

R(6) 

<3 

R(4) 

84 (10) 

54 (17) 

R(16) 

< 5 

R(5) 

6 (2) 

R(3) 

67 (9) 

55 (18) 

R(15) 

5(2) 

R(4) 

6 (2) 

R(2) 

73 (10) 

64 (20) 

R(14) 


R(3) 

< 5 

R(l) 

70 (9) 

39 (12) 

R(13) 

5(2) 

R(2) 


R(0) 

68 (9) 


R(12) 

< 4 

R(l) 

■ * • 

P(l) 

65 (9) 

11 (4) 

R(ll) 

6(2) 

R(0) 

< 6 

P(2) 

72 (10) 

30 (10) 

R(io) 

6(2) 

P(l) 

< 3 

P(3) 

80 (11) 

21 (7) 

R(9) 

7(2) 

P(2) 

< 3 

P(4) 

69 (9) 

82 (26) 

R(8) 


’ ' ’ 


P(5) 

84 (11) 

73 (23) 

R(7) 

12 (3) 



P(6) 

85 (10) 

65 (20) 

R(6) 

13 (2) 



P(7) 

85 (11) 

90 (28) 

R(5) 

19 (2) 



P(8) 

85 (11) 

71 (22) 

R(4) 

17 (2) 



P(9) 

85 (11) 

69 (21) 

R(3) 

21 (2) 



P(10) 

86 (10) 

84 (26) 

R(2) 

18(2) 



P(ll) 

79 (10) 

61 (18) 

R(i) 




P(12) 

79 (11) 

56 (17) 

R(0) 

24(2) 



P(13) 

88 (12) 

43 (13) 

P(i 

22 (2) 



P(14) 

69 (9) 

64 (19) 

P(2) 

25 (2) 



P(15) 

72 (10) 

53 (16) 

P(3) 

18 (2) 



P(16) 



P(4) 

15(2) 




Table 2 

Physical Parameters From Curve of Growth and Rotation Diagram 


Trot 


O 

o 

Vlsr 

Origin 

K 

km s" 1 

10 18 cm 2 

km s“ 1 


19^ 

15±l 

13+f 

1.6 

6 

5±3 

circumstellar disk+foreground; from l3 CO 

1.3 

10 

5±3 

as above, but alternative bo 

1.0 

19 

5±3 

as above, but alternative bo 

nil 

0.7 

96 

5±3 

as above, but alternative bo 

300l?“ 

250^f™ 

1.3 

1.0/0. 7 

2.4 

2.8 

5±3 

5±3 

disk; from 13 CO 

as above, but alternative bo 

250lf°° 

< 32 

O^OiO.OS 6 

28±6 

12 CO red wing 

500±250 

< 12 

~ 0.06 6 

-5/15 

blue and red wings resolved 12 CO main 


a assuming N( 12 CO)/N( 13 CO)=80 
assuming optically thin absorption 


nificant amount of warm gas (T ~ 250 K) are present 
within ~ 3 km s“* of the systemic velocity (Table 2). The 
12 CO lines show that warm gas at T ~ 250 K is also 
present at highly red shifted velocities (20 — 100 km s” 1 ), 
but at a factor 10 lower column. A small amount of warm 
gas is present at low red and blue shifted velocities as well 
(within 10 km s^ 1 ). As further described in §4.1, the gas 
components at the red shifted and systemic velocities can 
be fitted within the same physical model of a contracting 
disk, but the origin of the warm gas at low blue shifted 
velocities is more difficult to explain. 

3.2. Solid CO 

The broad absorption feature between 2122-2149 cm -1 
(Fig. 1) can be entirely attributed to the stretching vibra- 


tion mode of 12 CO in circumstellar ices (§4.3). The high 
spectral resolution allows, for the first time, to unambigu- 
ously separate the gas phase CO lines from the solid state 
absorption and study the solid CO band profile in great 
detail. In accordance with previous, low resolution stud- 
ies (Chiar et al. 1998; Teixeira et al. 1998), a distinct 
narrow feature is seen at 2140 cm -1 , and a significantly 
broader component at longer wavelengths. Our data how- 
ever, indicates the presence of a new, third component on 
the blue side, separate from the narrow 2140 cm -1 feature, 
most notable by a change of the blue slope at 2142 cm -1 
(Fig. 1). 

In order to explain the shape of this CO absorption pro- 
file, we have taken laboratory experiments of solid CO 
from the literature (Sandford et al. 1988; Schmitt, Green- 
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berg, & Grim 1989; Palumbo & Strazzulla 1993; Trotta 
1996; Ehrenfreund et al. 1997; Elsila, Allamandola, & 
Sandford 1997; Baratta & Palumbo 1998; Teixeira et al. 
1998). We determined the peak position and width of 
the laboratory profiles as a function of ice composition, 
temperature, and cosmic ray bombardment intensity. Ad- 
ditionally, for CO-rich ices we used optical constants to 
calculate the absorption profile as a function of particle 
shape in the small particle, Rayleigh limit (for details, see 
Ehrenfreund et al. 1997). The absorption profile of C0- 
poor ices (concentration <30%) is not affected by these 
particle shape effects. Changing these various parameters 
gives a wide variety of peak positions (2135-2144 cm" 1 ) 
and widths (FWHM=1.5-14 cm -1 ) in the laboratory. We 
limit ourselves here by identifying general trends, as sum- 
marized in Fig. 4. These trends are then related to the 
Gaussian peak position and width of the three aforemen- 
tioned components observed toward LI 489 IRS (also given 
in Fig. 4). 



Fig. 4. — Diagram showing schematically the effect of several as- 
trophysically relevant parameters on the peak position and width 
of the solid 12 CO stretching mode, as determined from laboratory 
experiments. Ellipses indicate ices of various composition. For CO^ 
rich ices, the effect of grain shapes is indicated by arrows for spheres 
and ellipsoids (‘CDE’). For polar ices, the effect (magnitude and di- 
rection) of heating is given by the arrow. ‘CH3OH+ CIT means a 
CH3OH or a CH30H:H20 ice irradiated by energetic particles sim- 
ulating cosmic rays. The dots with error bars represent the observed 
peak position and width of the absorption components detected to- 
ward L1489 IRS. 

The absorption band of a thin film of pure, solid CO at 
T =10 K peaks at ~2139 cm -1 , and is very narrow (~2 
cm" 1 ). This clearly does not fit any of the components 
observed toward L1489 IRS. However, for ellipsoidally 
shaped grains (more precisely, a distribution of ellipsoidal 
shapes, ‘CDE’; Bohren & Huffman 1983), this sensitive 
strong band shifts to shorter wavelengths and becomes 
broader. Now it exactly fits the 2140 cm _1 component 
observed toward L1489 IRS, both in peak position and 
width (Fig. 4). In view of a recent controversy on optical 
constants (Ehrenfreund et al. 1997), it is worth noting 
that good fits are obtained with optical constants from 
the works of Ehrenfreund et al. (1997) and Baratta & 
Palumbo (1998). The optical constants of Trotta (1996) 


do not induce strong particle shape effects and therefore 
pure CO does not provide a good fit to L1489 IRS for any 
particle shape. 

Broadening of the laboratory profile, in order to fit the 
2140 cm" 1 feature in L1489 IRS, is also achieved by adding 
a small amount of C0 2 , O2, or H 2 0 molecules. To avoid 
a too large broadening, and minimize aforementioned par- 
ticle shape effects, this mixture needs be diluted in N 2 . 
This astrophysically relevant molecule does not broaden 
the feature, and gives a small blue shift (Ehrenfreund et 
al. 1997; Elsila et al. 1997), required to fit the 2140 cm" 1 
feature in L1489 IRS. Thus, both this mixture, as well as 
ellipsoidally shaped pure CO ice grains, provide good fits 
to the central 2140 cm -1 feature. 

Now, the width significantly increases, and the peak 
shifts to longer wavelengths by diluting CO in a mixture 
of molecules with large dipole moments such as H2O or 
CH3OH (Sandford et al. 1988; Tielens et al. 1991). This 
particular behavior is needed to fit the broad long wave- 
length wing seen toward L1489 IRS. Solid CH3OH has an 
abundance less than a few percent of solid H 2 0 toward 
low mass objects (Chiar, Adamson, & Whittet 1996). 
H2O seems the best dilutant, because of its large inter- 
stellar abundance. The use of H20:CO mixtures requires 
interesting additional constraints. Although a low temper- 
ature, unprocessed H2 O-rich ice does provide a good fit to 
the red wing, it can be excluded based on the presence of 
a prominent second absorption at ~2150 cm" 1 in the lab- 
oratory, which is clearly absent toward L1489 IRS. This 
second peak is caused by CO molecules located in pock- 
ets in an amorphous ice. These CO molecules are weakly 
bound, and the ~2150 cm" 1 peak disappears rapidly at 
higher T or as a result of cosmic ray hits (Sandford et al. 
1988). Thus, the H 2 0 ice responsible for the long wave- 
length wing toward LI 489 IRS must be thermally ( T > 50 
K) or energetically processed. 

The blue wing seen at ~2143 cm" 1 toward L1489 IRS 
can only be explained by an apolar ice. Adding a signifi- 
cant amount of C0 2 to a CO ice (C0 2 /CO> 1) results in 
the blue shift and broadening required to fit the observed 
wing. Somewhat less C0 2 is needed (C0 2 /CO~ 0.5) when 
a large amount of 0 2 is present. N 2 may be added as 
well, but is not essential except as a dilutant to reduce 
the effects of particle shape. A good fit is obtained by the 
mixture N 2 :O2:CO 2 :CO=l:5:0.5:l, as proposed in Elsila et 
al. (1997). If C0 2 /CO< 0.5 the band peaks at too high 
wavelength. In CO-rich ices this problem can however be 
overcome by particle shape effects (Fig. 4). In view of this 
effect it is not possible to constrain the relative molecular 
abundances of this interstellar component in more detail, 
but it is clear that an apolax C0 2 or 0 2 ice is needed, 
different from the distinct 2140 cm” 1 feature. 

A three component fit to the entire CO ice band of L1489 
IRS is shown in Fig. 5. Although this is not a unique fit, 
it does obey the global trends that we identified in the 
laboratory experiments. 

Finally, the solid CO column density is derived by di- 
viding the integrated optical depth over the band strength 
A. We take A = 1.1 x 10" 17 cm molecule” 1 independent 
of ice composition (Gerakines et al. 1995), and thus find 
/'/(solid CO)=6.5xl0 17 cm” 2 . The matin source of uncer- 
tainty here is in A, which is about 10%. CO in polar ices 
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contributes 3.5 x 10 17 cm 2 to the total column and the ap- 
olar components at 2140 and 2143 cm -1 contribute each 
1.5xl0 17 cm -2 . 



2150 2140 2130 2120 

v (cm -1 ) 


Fig. 5- — Observed CO absorption band of L1489 IRS, fitted 
by a combination of three laboratory ices: a polar H 2 0:CO=4:l 
(T=50 K; dotted) to account for the long wavelength wing, an apo- 
lar N 2 :O 2 :CO 2 :CO=l:5:0.5:l (T=10 K; dashed) for the short wave- 
length wing, and a pure CO ice (’CDE’ shape; T=10 K; solid) that 
fits the central peak. The thick, smooth gray line is the sum of these 
components. 


4. DISCUSSION 
4.1. An Infalling Disk 

The astrophysical meaning of the apparent two compo- 
nent temperature structure seen in the 13 CO rotation dia- 
gram (Fig. 3) requires further investigation. For high mass 
protostars it was found that similar rotation diagrams can 
be ‘mimicked’ by power law models of spherical envelopes 
(van der Tak et al. 2000). 

For L1489 IRS, the detection of molecular gas at a range 
of temperatures and red shifted velocities could indicate 
the presence of infalling gas at a range of radii from the 
protostar. Indeed, a 2000 AU radius contracting, disk- 
like structure was found in millimeter wave interferometer 
data (Hogerheijde & Sandell 2000). In a detailed follow- 
up study, Hogerheijde (2001) adopts a flared-disk model 
based on Chiang & Goldreich (1997) with a radial power- 
law distribution for the temperature 

T = 34(i?/1000 AU) -0 ’ 4 K, (1) 

and a density distribution that has a power-law drop-off 
with radius and a vertical exponential drop-off with scale 
height h 

p(R, z) = po(R/l 000 AU) - 1 - 5 e -z2 / ft2 kg cm- 3 . (2) 

The scale-height h is assumed to be a simple function of 
R } h = Rf 2. An inward-directed radial velocity field de- 
scribed as 

v; n = 1.3(f?/100 AU) -0,5 km s" 1 (3) 

is inferred, in addition to Keplerian rotation around a 
0.65 M© central star. 

Can this contracting disk model, based on (sub-) mil- 
limeter emission observations with angular resolution of 


4-8", reproduce the observed infrared CO absorption line 
profiles measured along a pencil beam? The absorption 
lines are modeled with the radiative transfer code of Hoger- 
heijde & van der Tak (2000); the high densities in the 
disk ensure LTE excitation for the lines involved, and line 
trapping is neglected in the excitation calculation. The 
model spectra include dust opacity at a standard gas/ dust 
ratio, as well as a N{ CO) = 1 x 10 18 cm -2 column of 
cold foreground material (15 K; Hogerheijde 2001); both 
factors do not affect the spectra in any significant way. 
The calculated spectrum is convolved with a Gaussian of 
FWHM=12 km s" 1 , which is the NIRSPEC instrumental 
resolution. 

We find that, while keeping all other parameters the 
same as in Hogerheijde (2001), the assumed density pro- 
file sensitively influences the wings of the 12 CO lines. This 
is enhanced by the fact that we are observing the flared 
disk of L1489 IRS at an inclination between 60° and < 90° 
(cf., Padgett et al. 1999), and the pencil beam crosses the 
disk at a few scale heights. Small changes in the density 
profile, for example induced by the thermal structure, have 
a large effect on the absorption line profile. In the model of 
Hogerheijde (2001) the scale height increases linearly with 
distance from the star, and thus the density p{V) along the 
line of sight l follows the density in the mid-plane (Eq. 2) 
reduced by a factor e - 4 / t&n2 (°), with a being the inclina- 
tion. Here, we include the effect of density variations, or 
deviations from the adopted scale height h = R/ 2, by re- 
laxing the values of the density along the line of sight l by 
fitting p(l) = po(//1000 AU)” P to the data. 

This initial model successfully fits the peak velocity and 
depth of both high and low J 13 CO lines (Fig. 6). Its ro- 
tation diagram is quite different from that of the curve 
of growth analysis (Fig. 3), showing that rotation dia- 
grams must be interpreted with great care. Our model 
also matches the range of velocities observed in the red 
wings of the 12 CO lines, when taking p = 0.55 ± 0.15. 
This is a much shallower density profile compared to that 
derived from millimeter wave data (p = 1.5; Eq. 2) and 
indicates that the scale height increases more than lin- 
ear, i.e. the disk flares more than assumed in Hogerheijde 
(2001). With this result, it is possible to determine the 
important relation of disk scale height h(R) = a,R b as a 
function of R, but only if the disk inclination is a priori 
known. Unfortunately the inclination is not better con- 
strained than within the range of 60° and < 90° imposed 
by near-infrared data (Padgett et al. 1999). We can there- 
fore not distinguish between low and high values of a and 
corresponding high and low inclinations respectively. In 
either case, the total 12 CO column along the pencil beam 
is 1.2 x 10 19 cm -2 , with 58% of the CO mass at a tem- 
perature of T = 20 — 60 K, 15% at 60 — 90 K, and 27% 
at 60 — 90 K. This result is of importance in §4.3 in the 
interpretation of the solid CO observations, and in par- 
ticular in assessing the thermal history of ices. The total 
column of our model is in good agreement with the col- 
umn derived from the visual extinction (1.4 x 10 19 cm" 2 ; 
§3.1). It is also of the same order of magnitude as the 
total column through the mid-plane, calculated from dust 
and line emission (iV[CO] = 6 x 10 18 cm -2 ; Hogerheijde 
2001), and confirms the relatively edge-on orientation of 
the disk. 
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Wavenumber (cm -1 ) 

Fig. 6. — Comparison of the observed spectrum of L1489 IRS with an infalling disk model. The histograms in panels (a) and (b) represent 
the average of the observed 12 C0 P(6)-P(15) lines and all 13 CO lines respectively, corrected for the source and earth velocity (43 km s -1 ). 
The smooth thick gray line in these two panels is the collapsing disk model, averaged over the same lines. In panels (c) and (d) the same 
averaged observed spectra and model are plotted, but now assuming that in the model 1 % of the unextincted stellar continuum emission does 
not pass through the disk. Panel (e) shows the latter model (thick gray line) compared to a portion of the non-averaged L1489 IRS spectrum. 


However, apart from these successes, the 12 CO lines 
show that our infalling disk model produces too much 
warm gas at high velocities (Fig. 6). The 12 CO lines are a 
factor of 2.5 deeper, and, in contrast to 13 CO, they peak at 
a too high velocity (+10 km s" 1 ) with respect to the obser- 
vations. In principle, one could make the 12 C0 lines less 
deep by assuming that ~ 1% of the original, unextincted 
continuum flux (corresponding to 30% of the extincted 
continuum) reaches the slit without passing through the 
disk, by scattering on large grains. The shift in peak ve- 
locity however requires a solution of a more fundamental 
origin. Perhaps the infall velocity function is shallower, 
and the disk is more rotationally supported at lower radii. 
The amount of warm gas at high velocities can also be low- 
ered by assuming that only part of the disk participates in 
the high velocity inflow, such as a thin hot surface layer, or 
gas accelerated in magnetic field tubes directed from the 
inner disk to the stellar photosphere. Such a two compo- 
nent model is consistent with the rotation diagram derived 
from the curve of growth (Fig. 3), and also with the low 
observed mass accretion rate. If we take the inflow at face 
value, and assume that the entire disk participates, the 
mass accretion rate would be 1Q“ 6 M©, generating 7 Lq 
in accretion luminosity. The star’s Lboi is estimated at 3.7 
L© which also contains the stellar luminosity. It is there- 
fore indeed likely that the mass accretion onto the star is 
significantly lower, as is also traced through the lack of 


the hydrogen Pf/? emission line in our spectrum (2148.8 
cm" 1 ; Fig. 1) and the weakness of Br 7 emission (Muze- 
rolle, Hartmann, & Calvet 1998). 

4.2. Binarity? § 

An entirely different explanation for the line profiles may : 
lie in the possibility that L1489TRS is a protobinary sys- - 
tem. A protobinary nature of L1489 IRS is suggested by 
various pieces of evidence (Lucas, Blundell, & Roche 2000; 3 

Wood et al. 2001 and references therein). The presence 
of a quadrupolar outflow system is inferred from K band = 
polarization images, C 18 0 emission line profiles, Herbig- I 
Haro knots that are scattered throughout the L1489 IRS f 
environment, and a very complex near infrared scattered 
light pattern. Three dimensional models, in which the ax- 
isymmetry of the infalling circumstellar envelope is broken 
by multiple outflow cavities that are perpendicular to each 
other, are able to account for the observed morphology. 

The putative binary itself, however, has not been re- 
solved so far. An upper limit on the projected separation 
has been set at < 20 AU from near infrared images (Pad- 
gett et al. 1999). If the CO absorption line profile is 
in any way related to a binary system, then the large ob- 
served velocities (~23 km s" 1 ; §3.1) may indeed favor a 
close binary system. The line profile is then expected to 
vary on a time scale of a few months, which can easily be 
tested. In this case, much of the observed warm gets might 


in mill 
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be present in two small circumstellar disks, which are in a 
close orbit around each other. Some of the warm gas may 
also be present at low density in the central cavity created 
by the binary. The large column of cold gas may originate 
in the circumbinary disk. The binary extracts momentum 
from the 2000 AU circumbinary disk, setting up the in- 
ward motion seen in millimeter wave emission lines. We 
leave further investigation of this topic for future studies. 

4.3. The Origin and Evolution of Ices 

In order to establish if the solid CO observed toward 
L1489 IRS originates in foreground clouds or in a circum- 
stellar (or binary) disk, it is worth to compare with ices 
observed in lines of sight not affected by star formation. 
Observations of field stars obscured by intervening quies- 
cent material of the Taurus Molecular Cloud have revealed 
that solid CO is not present when the extinction A\ < 5 
(e.g. Teixeira & Emerson 1999). The solid CO toward 
L1489 IRS can therefore not be associated with foreground 
clouds, which have a gas column of 7V( 12 CO)= lxl0 ls 
cm" 2 (Hogerheijde 2001), corresponding to A v ~ 2. 

Thus, the solid CO must be present in the disk of L1489 
IRS. The absorption profile is intriguingly different from 
that seen in quiescent clouds. The broad red wing has 
a depth of ~30% with respect to the narrow 2140 cm -1 
peak, which is significantly more them toward all measured 
background stars (10%; Chiar et al. 1995). This may 
well be an effect of thermal processing along the L1489 
IRS line of sight, because the sublimation temperatures 
of polar and apolar ices, causing the broad and narrow 
features respectively, are very different (90 versus 18 K). 
However, a chemical origin of an increased abundance of 
polar ices in disks cannot be excluded, because the ap- 
parently edge-on system Elias 18 in the Taurus Molec- 
ular Cloud has an extremely large CO depletion factor 
(solid / [gas-bsolid] ^ 1 00% versus 7% for L1489 IRS), but a 
deep red ‘polar’ CO wing is present as well (Shuping et al. 
2001; Chiar et al. 1998). On the other hand, energetic 
processing may take place even in the cold disk of Elias 
18 (Whittet et al. 2001). Clearly, it is necessary to ob- 
servationally characterize the ices in circumstellar disks in 
much more detail. 

If for now we assume the sublimation scenario, we can 
do some general extrapolations which can be compared 
with the results of our gas phase study (§4.1). By scal- 
ing the long wavelength wing of solid CO of background 
field stars to that of L1489 IRS, we find that a column 
of 6xl0 17 cm" 2 of CO has evaporated from the apolar 
ice component in the part of the L1489 IRS disk along 
the pencil absorption beam where T < 90 K (the subli- 
mation temperature of polar ices). Then the column of 
solid CO that went from the quiescent cloud into building 
this part of the disk is 12.5 xlO 17 cm" 2 . Extrapolating 
this further, we use the observed CO depletion factor of 
30% toward field stars behind the Taurus Molecular Cloud 
(Chiar et al. 1995) to calculate that the original quiescent 
gas column must have been of the order of 3xl0 18 cm" 2 . 
Adding the evaporated column, the expected present day 
gas column at T < 90 K is 3.6xl0 18 cm -2 . This is of the 
same order of magnitude as the CO column below 90 K in 
our collapsing disk model (i\r[CO]=8.7 x 10 18 cm -2 ), which 
may indicate that no chemical change in the apolar/polar 


CO ice ratio and no significant additional depletion has 
occurred in the evolution from quiescent Taurus Molec- 
ular Cloud material to the formation of the LI 489 IRS 
disk. This contrasts strongly with the very large deple- 
tions found in the (older) disks of T Tauri stars (Dutrey 
et al. 1997). The low CO depletion along the pencil beam 
toward L1489 IRS (7%) and the supposed signs of thermal 
processing (see below) may be due to the fact that our line 
of sight does not cross the disk mid-plane, i.e. the system 
is not exactly edge-on. The ice processing we see takes 
place higher in the disk atmosphere, perhaps in the warm 
layer below the super-heated dust layer responsible for mil- 
limeter wave line emission (van Zadelhoff et al. 2001). It 
must be noted that in the model of Hogerheijde (2001) 
the gas temperatures are larger than 25 K, prohibiting the 
formation of apolar ices and large CO depletions anywhere 
in the disk. The observed presence of apolar CO ices thus 
indicates that, as already suggested in §4.1, the line of 
sight may cross the cold, rotationally supported disk in- 
terior not traced in the observations and infall model of 
Hogerheijde (2001). 

Apart from evaporation of apolar ices, other hints of 
thermal processing include the aforementioned absence of 
the 2150 cm" 1 absorption (§3.2), which occurs in cold un- 
processed polar CO ices but disappears at temperatures 
T > 50 K. Also, the blue apolar wing may be a conse- 
quence of thermal processing. If the central 2140 cm" 1 
peak is due to a mixture of 02? N2 and CO instead of 
pure CO (spectroscopically these cannot be distinguished), 
thermal or energetic processing (UV radiation from the 
ISRF, UV induced by H 2 cosmic ray collisions, or direct 
hits of cosmic rays) could efficiently produce C0 2 . This 
could cause the band to broaden and shift to the position 
of the observed blue wing. Chemical models indicate that 
energetic processing of molecules in disks takes place on 
a time scale of 10 6 yrs (Aikawa et al. 1999), which is 
somewhat longer than the age of the disk of LI 489 IRS 
(~ 5 x 10 5 yrs). This however applies to the disk mid- 
plane, and the time scale may well be shorter in the lower 
density higher disk layers that our observations of L1489 
IRS trace. A possible problem with the energetic process- 
ing interpretation is the absence of a feature adjacent to 
the short wavelength side of the CO ice band, usually at- 
tributed to energetically produced C=N bondings (Whit- 
tet et al. 2001). Another spectroscopic tracer of thermal 
processing is the signature of crystallization in the band 
profiles of H 2 0 and C0 2 ices. Our infalling disk model 
predicts that only 15% of the gas is within the tempera- 
ture range at which ices crystallize (60-90 K), and thus 
crystallization is not expected to play a significant role in 
the disk of L1489 IRS. This model prediction can be tested 
with future high quality H 2 0 and C0 2 spectra of L1489 

ms. 

In summary, several pieces of evidence indicate that the 
CO ices in the disk of L1489 IRS have experienced thermal 
or energetic processing. The strongest arguments are the 
low depletion factor and the low ratio of apolar to polar 
ices with respect to the quiescent Taurus Molecular Cloud 
material. This may be explained by the fact that the disk 
of L1489 IRS is seen under an angle, and our pencil ab- 
sorption beam traces the warm upper disk layers. 
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5. SUMMARY AND FUTURE WORK 

We have shown that valuable and unique information is 
obtained from high resolution spectroscopy of the CO fun- 
damental at 4.7 fim toward the low mass class I protostar 
L1489 IRS in the Taurus Molecular Cloud. At a resolution 
of R = 25,000 (12 km s” 1 ) this object shows a multitude 
of deep ro-vibrational absorption lines of 12 CO, as well as 
13 CO and C 18 0. The isotopes trace large columns of warm 
and cold gas in the circumstellar disk at or within 3 km s” 1 
of the systemic velocity, while the 12 CO line profiles show 
warm gas that is red shifted at a range of velocities of 
up to 100 km s" 1 . Both the line depth of the isotopes 
and the extent of the red shifted warm gas seen in 12 CO 
are well explained by an infalling flared disk model with 
power laws for the temperature, infall velocity, and density 
(a small column of blue shifted gas seen in the 12 CO line 
wing however remains unexplained). These observations 
show that the inward motions inferred on scales of several 
hundred AU through millimeter wave interferometry con- 
tinue to within 0.1 AU of the star, where the velocity model 
of Hogerheijde (2001) predicts inward velocities exceed- 
ing several tens of km s" 1 . A detailed comparison of our 
power law infall model however overestimates the amount 
of warm, high velocity infalling gas. Much of this gas must 
therefore be rotationally supported, and only a thin disk 
surface layer is infalling, or gas is accelerated along mag- 
netic fields in the inner parts of the disk. High spatial 
resolution millimeter wave observations (with ALMA) are 
needed to test our model, e.g. to refine the determination 
of the velocity field in the inner disk parts and the de- 
pendence of the disk scale height on radius h(R). Finally, 
high spatial resolution infrared interferometer observations 
would be able to see if LI 489 IRS is a close binary system 
(< 20 AU), which is essential to assess the importance of 
this aspect to the observed line profiles and disk evolution. 

In the same spectrum, a deep CO ice band is seen toward 
L1489 IRS, which we conclude to originate in the disk. The 
high signal-to-noise and spectral resolution allowed us to 
separate it from circumstellar CO absorption lines, and 
to do a detailed comparison with laboratory mixtures. It 
is found, for the first time, that besides the well known 


long wavelength wing of polar H 2 0:CO ices and the cen- 
tral narrow absorption peak at 2140 cm" 1 , an additional, 
separate component is present on the blue side of the ice 
band. Both the central peak and the blue shifted compo- 
nent are due to CO in apolar ices, one due to pure CO and 
one due to a mixture of C0 2 and perhaps 0 2 and N 2 . The 
relatively large depth of the long wavelength, polar wing 
relative to the apolar components may suggest that the 
ices in the disk of LI 489 IRS are thermally processed with 
respect to the quiescent Taurus Molecular Cloud material, 
which is supported by the low depletion factor (~7%) in 
this line of sight. This is likely due to the fact that we see 
the disk under an angle, and our pencil absorption beam 
traces the upper disk layers. 

The present work shows that high spectral resolution 4.7 
lim observations are a great tool to better understand pro- 
tostellar disks, which define the initial conditions of planet 
and comet formation, both in the solid and gas phase. This 
is an exploratory study and needs to be followed up by ob- j 
serving a larger sample of protostellar disks at 4.7 fim to j 
investigate the influence of parameters such as foreground 
contribution, disk inclination, age, and binarity on gas and 
ice band profiles. 
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